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ABSTRACT 


Traction  tests  were  performed  on  RP1,  a  common  kerosene  based  rocket 
propellant.  Traction  data  on  this  fluid  are  required  for  purposes  of 
turbopump  bearing  design,  using  codes  such  as  SHABERTH.  To  obtain  the 
traction  data,  an  existing  twin  disc  machine  was  used,  operating  under  the 
side  slip  mode  and  using  elliptical  contacts.  The  range  of  test  variables 
were:  contact  peak  Hertz  stress  from  1.0  to  2.0  GPa,  disc  surface  speed  from 
10  to  50  m/s,  fluid  inlet  temperature  from  30  to  70  °C,  and  with  a  contact 
aspect  ratio  of  1.7. 

The  resulting  traction  curves  were  reduced  to  fundamental  fluid  property 
parameters  using  the  Johnson  and  Tevaarwerk  traction  model.  Theoretical 
traction  predictions  were  performed  by  back  substitution  of  the  fundamental 
properties  into  the  traction  model.  Comparison  of  the  predicted  with  the 
measured  curves  gives  a  high  degree  of  confidence  in  the  correctness  of  the 
traction  model.  For  purposes  of  input  to  the  NASA  SHABERTH  program,  the 
traction  model  was  next  used  to  predict  the  expected  traction  of  RP1  under 
line  contact  conditions. 
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1-0  INTRODUCTION 


Traction  or  friction  plays  a  major  role  in  today's  technological  society 
in  that  it  holds  one  of  the  keys  to  reduce  our  overall  energy  consumption,  and 
thereby  the  dependence  on  unreliable  sources  of  this  energy.  Friction  and 
traction  indicate  the  resistance  to  relative  motion  of  two  'contacting' 
bodies.  The  term  traction  and  friction  have  the  same  meaning  in  a 
tribological  sense,  however  friction  is  used  when  this  resistance  is 
undesirable  and  traction  is  used  when  it  is  desirable.  The  mechanical 
components  in  which  friction  and  traction  are  important  are  rolling  element 
bearings,  gears,  cam  and  tappets,  and  traction  drives.  Because  these  devices 
almost  always  operate  in  a  wet  or  fluid  lubricated  environment,  the  traction 
or  friction  is  mostly  governed  by  the  particular  fluid  that  is  used.  In  the 
first  three  devices  friction  is  the  key  source  of  inefficiency  and  because  of 
the  multitude  of  bearings  and  gears  in  service,  a  small  reduction  in  these 
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losses  can  amount  to  phenomenal  savings  in  energy.  Rolling  element  bearings 
actually  have  rather  interesting  requirements  for  friction  or  traction  in  that 
at  low  to  medium  speeds  friction  should  be  low,  but  at  high  speeds  traction 
should  be  high  to  ensure  that  the  rolling  elements  operate  at  the  correct 
velocities. 

In  most  mechanical  devices  where  friction  or  traction  plays  an  important 
role,  the  lubricating  fluid  is  chosen  to  optimize  the  energy  efficiency  of  the 
device,  however,  this  is  not  possible  in  all  cases.  For  example,  the  bearings 
used  in  rocket  motor  turbopumps  are  lubricated  with  the  process  liquid,  may 
this  be  a  mineral-based  low-viscosity  fuel  such  as  RP1,  or  a  cryogenic 
substance  such  as  liquid  oxygen  or  hydrogen.  For  applications  such  as  these 
it  is  critically  important  that  the  rheological  properties  of  the  working 
flui^  be  known  and  modeled  such  that  a  good  deal  of  confidence  exists  about 
using  the  data  beyond  the  range  of  observation.  The  data  thus  obtained  can 
then  be  used  in  bearing  design  codes  such  as  SHABERTH  to  design  for  adequate 
bearing  performance  for  a  given  mission. 

1-1  Prior  Traction  Investigations 

There  has  been  a  lot  of  activity  in  the  area  of  traction  research,  both 
in  the  past  and  recently.  Notable  contributions  have  come  from  Clark  et  al 
[1951],  Smith  [1965],  bmith  et  al  [1973],  Johnson  and  Cameron  [1967],  Niemann 
and  Stoessel  [1971],  end  more  recently  Johnson  and  Roberts  [1974],  and  Johnson 
and  Tevaarwerk  [1977].  Some  of  these  investigations  were  strictly 
experimental  in  nature  and  aimed  it  obtaining  traction  design  data,  while 
others  were  aimed  at  understanding  the  traction  phenomena  so  that  rheological 
models  could  be  formulated.  This  latter  research  is  of  course  ultimately 
aimed  at  relating  fluid  molecular  properties  to  traction  properties.  Research 
by  Johnson  and  Tevaarwerk  [1977],  Daniels  [1978],  Hirst  and  Moore  [1980],  and 
Alsaad  et  al  [1978]  is  directed  more  towards  this  purpose.  The  reader  is 
referred  to  an  excellent  review  by  Johnson,  [1978]  for  further  aspects  of  this 
topic . 

Current  understanding  of  traction  has  led  to  traction  models  that 
describe  the  fluid  shear  behavior  in  terms  of  an  elastic  and  a  dissipative 
element.  For  purposes  of  mathematical  tractability  this  dissipative  element 
is  taken  to  be  plastic  like  in  nature.  This  gives  an  adequate  description  of 
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the  fluid  behavior  at  conditions  such  as  those  encountered  in  rolling  element 
bearings.  This  work  has  now  been  further  expanded  by  developing  a  simple 
method  to  correct  for  thermal  effects  due  to  spin,  Tevaarwerk  [1981],  and  an 
overall  thermal  traction  study,  Tevaarwerk  [1980]. 

As  with  all  models,  however,  their  usefulness  is  severely  restricted  if 
inadequate  input  data  is  available  to  the  designer.  This  is  especially  so  if 
fluids  are  used  that  are  not  normally  desirable  from  a  lubrication  point  of 
view,  but  are  selected  on  the  basis  of  other  features.  RP1  is  such  a  fluid. 
Being  a  rocket  propellent,  based  on  a  high  viscosity  kerosene,  it  does  not 
have  the  traction  data  base  that  other  more  lubricous  fluids  have.  Therefore, 
it  was  decided  to  test  its  traction  behavior  experimentally  and  to  use 
existing  rheological  models  to  predict  its  traction  under  conditions  different 
from  those  in  the  experiments. 

The  investigation  reported  herein  was  performed  by  Battelle  Memorial 
Institute,  of  Columbus,  Ohio,  under  contract  to  the  NASA  Lewis  Research  Center 
through  the  Metals  and  Ceramics  Information  Center,  a  DOD  "Information 
Analysis  Center",  Item  0001CJ,  Control  No.  DLA  900-83-C-1744.  The  NASA 
technical  project  manager  was  C.  Woods  of  the  Launch  Vehicle  Technology  Branch 
in  the  Space  Propulsion  Division. 
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2-0  EXPERIMENTS 

The  various  traction  experiments  were  carried  out  on  an  existing  twin 
disc  test  facility,  shown  in  Figures  2-1,  2-2  and  2-3.  The  test  facility  was 
also  used  for  the  traction  test  results  as  reported  in  Tevaarwerk  [1983]  and 
[1985].  Traction  curves  are  obtained  by  using  the  side  slip  technique  so 
that  large  traction  transfer  can  be  measured  without  the  need  for  a  large 
motor/generator  set.  It  also  has  the  advantage  that  the  effects  of  bearing 
friction  and  rolling  traction  can  be  corrected  for  by  the  simple  measurement 
of  the  longitudinal  slip.  This  is  particularly  important  when  the  expected 
levels  of  fluid  friction  are  low,  as  is  the  case  with  the  lubricant  tested 
here.  With  longitudinal  traction  measurements  the  correction  for  bearing 
friction  and  rolling  traction  is  often  of  the  same  magnitude  as  the  measured 
traction  with  these  lubricants. 

2-1  Description  of  Twin  Disc  Machine. 

For  an  extensive  description  of  twin  disc  traction  testers  the  reader  is 
referred  to  the  literature;  Smith  [1965]  and  Johnson  and  Roberts  [1974]. 
Basically  the  machine  consists  of  two  discs,  called  the  toroid  (upper)  and 
disc  (lower).  The  lower  disc  is  mounted  in  rolling  element  bearings  through 
shafts  and  the  only  degree  of  freedom  is  one  of  rotation  about  this  axis.  The 
lower  disc  always  has  a  transverse  radius  of  curvature  of  infinity.  The  upper 
disc  is  contained  in  bearings  that  are  mounted  in  the  upper  assembly.  This 
upper  assembly  is  suspended  with  elastic  hinges  such  that  only  direct  normal 
motion  and  axial  motion  is  possible.  The  assembly  will  however  always  stay 
horizontal.  The  toroid  (or  upper)  has  curvatures  to  produce  the  desired 
contact  geometry  is  arrived  at.  The  upper  assembly  is  constructed  such  that 
the  toroidal  axis  can  be  tilted  relative  to  the  horizontal  plane  so  as  to 
introduce  spin  on  the  contact.  It  can  also  be  skewed  about  the  normal  to  the 
contact  to  introduce  a  side  slip  velocity. 

In  order  to  achieve  the  range  of  pressures  required  for  the  traction 
data,  a  set  of  discs  with  a  nominal  aspect  ratio  of  1.7  is  used.  These  discs 
are  made  of  AISI-01  steel,  hardened  to  7.00  GPa,  ground  and  polished  to  a 
surface  finish  of  less  than  .05  p m  RMS  and  with  an  out-of-roundness  error  of 
less  than  5  ^m.  Between  tests  the  discs  are  inspected  for  surface  damage  and, 
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if  needed,  reground  and  polished  to  bring  them  back  up  to  specifications.  The 
required  normal  load,  obtained  by  a  dead  loading  technique,  can  be  calculated 
from  the  Hertz  theory  for  elastic  bodies  in  contact.  The  maximum  contact 
normal  stress  P0  is  given  by: 


(2-1) 


1  3  3  FZE  * 2 

Mx  K2 


where, 

Fz  =  contact  normal  load,  [N] 

P0  =  Hertzian  contact  stress,  [Pa] 

E'  =  composite  elastic  modulus,  [Pa] 

(231  GPa  for  steel) 

Re  =  equivalent  disc  radius,  [m] 

=  (1/RX  +  1  /Ry ) _1 ,  and 
x  =  Hertzian  contact  shape  factors.  [-] 


2-2  Instrumentation  of  the  Traction  Tester 

By  suitably  instrumenting  the  disc  machine  the  relevant  experimental 
parameters  can  be  measured.  In  the  experiments  reported  here  the  sideslip, 
sideways  traction  force,  disc  surface  temperature,  rolling  velocity,  and  the 
amount  of  longitudinal  slip  are  measured.  The  technique  of  measuring  each  of 
these  variables  will  be  discussed  next. 

2-2-1  The  Measurement  of  the  Traction  Force 

A  ring  dynamometer  type  load  cell  is  used  to  measure  the  side  slip  force 
of  the  upper  toroid  assembly.  The  electrical  signal  from  the  load  cell  is 
conditioned  for  noise  and  amplified  using  common  mode  rejection  techniques. 
The  gain  on  the  amplifier  is  adjusted  so  that  a  good  range  on  the  signal  is 
measured  for  each  test.  Calibration  of  the  load  cell  is  done  in  situ  by  dead 
loading.  This  calibration  is  checked  periodically. 
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2-2-2  The  Measurement  of  Side  Slip 


The  skew  angle  is  measured  by  using  a  direct  current  displacement 
transducer  on  the  upper  assembly  and  thereby  measuring  the  rotation  angle  of 
this  assembly.  This  skew  angle  gives  the  amount  of  side  slip/roll  ratio 
through  the  relationship: 


(2-2) 


Av/U  =  tan  (p)  , 


where, 

rad] 
m/s] 
m/s] 


p  =  side  slip  angle, 

Av  =  side  slip  velocity,  and 
U  =  rolling  velocity. 


By  measuring  the  amount  of  skew  with  the  displacement  transducer  the  side 
slip/roll  ratio  is  obtained  directly  through: 

(2-3)  av/U  =  k'Ax  . 

where  ax  is  the  displacement  of  the  transducer  and  K  a  scale  factor.  The 
electrical  output  of  the  displacement  transducer  is  filtered  in  an  R-C  network 
to  provide  a  low-pass  signal.  The  maximum  frequency  response  of  the  R-C 
network  is  0.1  s.  The  displacement  transducer  is  calibrated  by  rotating  the 
upper  assembly  through  a  known  angle  and  then  calculating  the  amount  of  side 
slip  for  this  angle. 

2-2-3  The  Measurement  of  Rolling  Velocity  and  Longitudinal  Slip 

The  rotational  velocities  of  both  the  disc  and  the  toroid  are  measured  by 
using  a  MC  6840  frequency  counter.  The  proximity  probes  for  the  control  of 
the  counters  are  mounted  near  the  support  bearings,  see  Fig.  2-3.  Each  shaft 
has  a  simple  protrusion  on  it  that  produces  a  pulse  once  per  revolution.  This 
pulse  is  used  to  turn  countdown  counters  on  when  controlled  to  do  so.  The 
next  successive  pulse  from  the  proximity  probe  stops  the  countdown.  When  each 
of  the  two  counters  has  completed  its  cycle,  a  flag  is  set  and  the  register 
contents  are  read.  From  the  difference  between  the  original  and  the  final 
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contents  of  the  registers,  the  amount  of  time  for  a  single  revolution  of  each 
shaft  is  calculated.  From  this  the  angular  speed  for  each  shaft  is 
determined,  and  hence  the  peripheral  velocity  of  each  roller  can  be  calculated 
if  the  roller  radius  is  known.  The  countdown  frequency  on  the  counters  is 
selected  such  that  an  accuracy  of  at  least  1:  10000  is  obtained.  By  employing 
this  method  the  slip  of  the  rollers  can  be  measured  every  five  or  six 
revolutions  of  the  shafts.  From  the  countdown  and  the  undeformed  roller  radii 
the  slip  is  calculated  as  follows: 

,o  _  o  U2  -  U1  _  r,(  R2  C2  w  (  R2  ,  C2  ^  r  n 

(2_4)  s  ■ 2  uttut  ■  2V  RT  “  rr )  /  v  m  +  rr  /  • 

This  can  be  approximated  if  the  amount  of  slip  is  small  by  the  following: 

(2-S)  s  =  (  R?  "  U?)1 


where, 


s  =  roller  slip, 

U  =  rolling  velocity, 

C  =  number  of  count  down  pulses,  and 
R  =  roller  radius. 


m/s] 

[-] 

m] 


Suffix  1  denotes  the  toroid  and  suffix  2  denotes  the  disc.  It  should  be 
stressed  that  the  actual  amount  of  slip  will  be  somewhat  different  when  two 
different  roller  materials  are  in  contact  because  of  the  deformation  of  the 
rollers  themselves.  The  errors  thus  introduced  can  be  removed  by  measuring 
the  rolling  slip  at  very  low  rolling  speeds. 


2-2-4  The  Measurement  of  the  Disc  Surface  Temperature 


In  the  analysis  and  reduction  of  the  test  data  it  is  important  that  the 
disc  temperature  be  used  as  the  reference  inlet  temperature  for  the  film 
thickness.  This  temperature  can  be  measured  by  embedding  a  thermocouple 
directly  below  the  surface  of  the  disc  and  then  to  take  this  signal  out 
through  slip  rings.  This  method  is  however  not  very  practical  when  a  number 
of  different  discs  are  involved  and  is  also  very  costly  from  an  installation 
point  of  view.  With  care  the  surface  temperature  can  also  be  measured  by 
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using  a  trailing  thermocouple  that  rides  on  the  disc  surface.  The 
disadvantage  of  this  technique  is  that  it  can  be  speed  sensitive  in  its 
response  because  of  frictional  heating. 

The  latter  technique  is  employed  here  and  care  is  taken  to  ensure  that 
the  contact  force  on  the  thermocouple  is  not  excessive.  A  temperature 
reference  bath  junction  ensures  that  the  same  reference  level  for  the 
thermocouple  is  used  at  all  times.  The  signal  from  the  thermocouple  is 
amplified  using  common  mode  rejection  techniques  to  minimize  the  influence  of 
electrical  noise  and  other  disturbances.  Calibration  is  done  by  the  boiling 
water  method  adjusted  for  sea  level  differences.  This  calibration  is  checked 
periodically.  Only  slight  deviations  are  encountered.  Because  of  the 
frictional  heating  at  the  junction/toroid  interface  a  variation  of  about  2  C° 
is  found  in  the  signal  between  stationary  discs  and  those  rotating  at  a 
surface  velocity  of  120  m/s.  The  overall  reproducibility  of  the  temperature 
measurement  is  better  than  2  C°. 

The  temperature  on  the  machine  is  regulated  through  the  use  of  heaters 
and  coolers  on  the  test  fluid.  This  test  fluid  is  allowed  to  circulate  freely 
before  the  start  of  a  test  series  in  order  to  bring  the  machine  up  to  a 
uniform  temperature.  No  specific  effort  is  made  however  to  maintain  a  given 
set  point  temperature  during  the  test  and  the  reason  for  this  will  become 
clear  in  the  analysis  of  the  results. 

2-3  Traction  Measurements 

Traction  curves  are  obtained  by  the  slow  rotation  of  the  upper  assembly 
from  a  positive  value  of  side  slip/roll  ratio  to  a  negative  value.  The 
signals  from  the  above  discussed  transducers  are  fed  into  a  digitizer  from 
where  they  are  led  into  a  data  logging  computer  for  plotting  and  storage  on 
magnetic  media  for  future  use.  The  computer  automatically  traces  the  force 
versus  slip  curve  on  the  screen.  By  reversing  the  direction  of  rotation  of 
the  machine,  a  duplicate  set  of  curves  can  be  obtained.  For  each  experiment 
500  data  points  are  taken  at  fixed  time  periods  of  0.2  seconds.  Multiple  data 
points  are  stored  as  separate  entries. 

After  the  completion  of  a  test  series  the  data  is  recalled  into  memory  of 
the  computer  and  further  manipulated.  This  manipulation  consists  of  the 
averaging  of  the  multiple  entries,  the  filling  in  of  any  gaps  in  the  data 
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through  forward  and  backward  interpolation,  the  comparison  of  the  traces  for 
the  forward  and  reverse  rolling  direction,  and  the  centering  of  the  traces 
about  the  center  lines.  After  the  centering  operation  the  data  is  smoothed  by 
a  * N '  point  averaging  technique  for  traction  points  after  the  peak  traction 
points.  For  storage  a  geometric  series  is  used  so  that  the  total  traction 
trace  is  now  represented  by  40  data  points  for  each  measured  variable.  These 
traces  are  stored  on  magnetic  media  and  used  for  further  manipulation  and  data 
extraction  at  a  later  point. 

2-4  Traction  Curves  For  RP1 


A  survey  of  the  literature  revealed  that  very  little  traction  data  is 
available  on  low  viscosity  fluids  such  as  RP1.  This  is  especially  so  in  the 
higher  speed  ranges  of  10  m/s  and  above.  Hence  it  was  decided  to  perform 
traction  tests  on  RP1  over  a  fairly  broad  speed,  pressure,  and  temperature 
range.  The  nominal  test  matrix  consisted  of  the  following: 


Contact  Pressure  P0: 

1,  1.3,  1.6,  and  2.0,  I 

GPa 

Rolling  speed  U: 

±  10,  ±  30,  ±50,  j 

'm/s' 

Temperature  60: 

25  and  75,  1 

!°c] 

Slip  AV/U: 

-0.05  to  0.20,  and  j 

[-] 

wl  ab 

Spin  -g-  : 

0.  | 

[-] 

This  test  matrix  leads  to  24  different  traction  test  curves.  The 
measured  value  being  the  traction  coefficient  Fy/Fz  as  a  function  of  slip. 
The  traction  data  measured  on  RP1  using  the  equipment  and  techniques  as 
described  above  are  shown  in  Appendix  I.  The  lettered  labels  in  each  curve 
are  used  only  for  identification  purposes  and  do  not  indicate  an  actual 
measured  data  point. 


FIGURE  2-1.  OVERVIEW  OF  HIGH  PERFORMANCE  SIDE  SLIP  TRACTION  TESTER 
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3-0  THEORETICAL  ANALYSIS 

In  order  to  understand  the  required  analysis  of  the  experimental  data  it 
will  be  helpful  to  consider  the  following  discussion  of  traction.  The  ability 
of  a  fluid  film,  trapped  under  high  pressure  in  the  elastically  deformed 
region  of  two  loaded  curved  elements,  to  transmit  a  tangential  force  from  one 
element  to  the  other  is  commonly  referred  to  as  friction  or  traction.  The 
magnitude  of  this  force  depends  on  several  variables  such  as  :1)  the  contact 
kinematic  conditions  of  slip,  spin  and  sideslip,  2)  the  fluid  present,  3) 
temperature,  pressure  and  operating  speeds.  We  will  examine  the  traction 
behavior  under  simple  slip  only. 

Under  conditions  of  increasing  slip  between  the  two  elements,  an 
increasing  traction  force  is  transmitted  up  to  a  certain  limit  at  which  point 
it  will  decrease  with  further  slip.  See  Figure  3-1. 


FIGURE  3-1.  TYPICAL  TRACTION  /SLIP  CURVE 

There  are  three  regions  identified  on  this  traction  curve  and  the 
behavior  in  each  of  these  regions  can  best  be  described  by  the  Deborah 
number  (De) .  For  a  simple  Maxwell  viscoelastic  model  this  number  is  the  ratio 
of  the  relaxation  time  and  the  mean  transit  time,  see  Johnson  and  Tevaarwerk, 
[1977]. 

(A)  The  linear  low  slip  region.  Thought  to  be  isothermal  in  nature,  it  is 
caused  by  the  shearing  of  a  linear  viscous  fluid  (low  De)  or  that  of  a  linear 
elastic  solid  (high  De) . 

(B)  The  nonlinear  region.  Still  isothermal  in  nature  but  now  the  viscous 
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element  responds  nonlinearly.  At  low  De  this  portion  of  the  traction  curve 
can  be  described  by  a  suitable  nonlinear  viscous  function  alone,  while  at  high 
De  a  linear  elastic  element  interacts  with  the  nonlinear  viscous  element. 

(C)  At  yet  higher  values  of  slip  the  traction  decreases  with  increasing 
slip  and  it  is  no  longer  possible  to  ignore  the  dissipative  shearing  and  the 
heat  that  it  generates  in  the  film.  Johnson  and  Cameron  [1967]  showed  that 
the  shear  plane  hypothesis  advanced  by  Smith  [1965]  does  account  for  most  of 
their  experimental  observations  in  this  region.  More  recently  Conry  et  al 
[1979]  have  shown  that  a  nonlinear  viscous  element  together  with  a  simple 
thermal  correction  can  also  describe  this  region. 


3-1  Isothermal  Traction  Analysis 


The  rheological  model  that  describes  the  traction  under  simple  slip  in 
all  three  regions  of  operation  fairly  well  is  the  J  &  T  traction  model  as 
presented  by  Johnson  and  Tevaarwerk  [1977]: 


(3-1) 


1 

G  dt 


+  F  (r) 


7  , 


where, 


r  =  shear  stress, 

G  =  shear  modulus, 

£  =  time,  and 
7=  shear  strain  rate. 


[Pa] 

[Pa] 

.s] 

[1/s] 


The  dissipative  function  F(r)  is  open  to  the  choice  of  the  researcher  to  fit 
the  observed  traction,  but  Johnson  and  Tevaarwerk  [1977]  found  that  the 
hyperbolic  sine: 


(3-2) 


l 


where, 

ts  =  non-linear  shear  stress  parameter.  [Pa], 
tj  =  local  viscosity  [Pa.s], 


described  all  of  their  experimental  results  in  regions  (A)  and  (B)  very  well. 
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At  higher  pressures  and  for  fluids  with  high  traction  coefficients  this 
dissipative  function  may  be  replaced  by  the  purely  plastic  behavior  of  the 
material : 

(3-3)  F (r)  =0  for  r  <  re  ;  F(r)  =  7  for  r  =  rc  , 


where, 

rc  =  limiting  shear  strength  of  the  fluid  [Pa]. 

Whether  the  perfectly  plastic  behavior  of  the  material  is  intrinsic  is  not 
clear.  Work  by  Johnson  and  Greenwood  [1980]  suggests  that  it  is  possibly  the 
result  of  thermal  behavior  of  the  sinh  model.  For  many  applications  the 
elastic/plastic  traction  model  is  adequate.  It  was  used  by  Tevaarwerk  and 
Johnson  [1979c]  and  Tevaarwerk  [1979b]  to  predict  traction  under  various 
conditions  of  slip  and  spin.  The  analysis  is  completely  isothermal  in  nature 
and  for  simple  slip  the  tract  ion  is  given  hy: 


(3-4) 


J4  -| 


tan^S  + 


(1  +  S2) 


where, 


S  = 


2  J1 

I  71 


J1  = 


lab 

IT 


AU 

T 


and  J4  = 


£ 


The  shear  strain  rate  in  the  fluid  was  taken  to  be  the  same  everywhere  in 
the  contact  and  assumed  to  be  constant  throughout  the  film  thickness.  Its 
magnitude  was  taken  to  be: 

o  aU 

7  =  T  * 

Equation  (3-4)  results  from  the  integration  of  stresses  caused  by  the 
shearing  of  an  elastic  element  of  pressure-independent  average-shear-modulus  G 
and  the  plastic  stresses  proportional  to  the  local  Hertzian  pressure.  The 
predicted  traction  from  an  elastic/plastic  model  compares  very  well  with  the 
experimentally  observed  values  for  combinations  of  slip  and  spin,  provided 
that  the  spin  or  slip  are  not  too  large.  Large  slip  or  spin  results  in  almost 
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purely  dissipative  stresses  over  the  contact  area  and  hence  non-isothermal 
behavior.  Traction  prediction  under  t^ese  conditions  is  still  possible  but 
the  thermal  effects  need  to  be  brought  into  the  picture.  Tevaarwerk  [1981a], 
[1979c]  presents  t  vo  techniques  for  calculating  such  spin  traction  curves. 

The  former  technique  requires  the  shape  of  the  traction  curve  in  the  large 
slip  regime  to  provide  a  simple  correction  to  the  isothermally  predicted  spin 
traction . 

3-2  Thermal  Traction  Analysis 

The  ability  to  separate  the  elastic  stresses  from  the  plastic  ones  can  be 
used  to  perform  a  thermal  traction  calculation.  The  analysis  presented  here 
follows  the  technique  outlined  by  Tevaarwerk  [1980]  and  [1983]. 

Equation  (3-4)  resulted  from  the  integration  of  isothermal  elastic  and 
plastic  stresses  over  the  contact  area  of  an  ellipse.  For  the  region  of 
contact  under  elastic  stress,  the  shear  energy  is  conserved  and  therefore  does 
not  give  rise  to  temperature  increases.  The  plastically  deforming  region, 
however,  is  non-conservative  and  a  temperature  rise  in  the  fluid  is  expected. 
This  will  lead  to  a  reduction  in  the  local  strength  of  the  fluid.  Equation 
(3-4)  may  therefore  be  better  written  in  its  elastic  and  plastic  parts: 


(3-5) 

J4 

=  +  J4e  - 

(3-6) 

J4e 

IS) 

* r 

n 

and 

t(1+S2)2 

(3-7) 

J4p 

=  —  [  tan'lS  + 

S(S2  -1) 

P 

JT 

Equation  (3-5)  may  now  be  modified  by  tc/t0  where  rc  is  the  average 
stress  under  thermal  conditions  and  r0  is  the  average  stress  under  isothermal 
conditions.  We  are  dealing  therefore  with  averaged  stresses  in  the 
plastic  region  of  the  contact  even  though  the  isothermal  stress  distribution 
is  according  to  the  Hertzian  pressure.  This  seemingly  contradictory 
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assumption  is  supported  by  theoretical  evidence  by  Tevaarwerk  [1979c].  The 
modified  equation  would  therefore  be: 


(3-8)  J4t  =  J4e  +  ^  J4p 


The  modification  term  tc/t0  can  be  found  from  a  thermal  balance  over 
the  contact  region  under  plastic  stress.  This  region  can  be  thought  of  as  a 
thermal  source  whose  heat  is  conducted/convected  away.  The  length  of  the 
source  is  a  function  of  the  location  of  the  onset  of  plastic  deformation  after 
the  initial  elastic  region,  however  in  this  simple  model  we  will  take  the 
source  length  to  be  "a"  where  this  is  the  semi  contact  length  in  the  running 
direction.  As  a  simple  thermal  balance  we  will  use  the  expression  reported  by 
Johnson  and  Cameron  [1967]  for  the  shear  plane  temperature: 


(3-9) 


where, 


F.  + 


f 


(3-10) 


[-] 


and, 


q  =  average  thermal  strength  of  the  source, 

0C  =  shearplane  temperature  in  the  contact, 

80  =  inlet  temperature  of  the  fluid, 

kf  =  thermal  conductivity  of  fluid  in  the  contact, 

h  =  central  film  thickness  in  the  contact, 

ks  =  thermal  conductivity  of  the  roller  material, 

Cs  =  specific  heat  of  the  roller  material, 

ps  =  density  of  the  roller  material,  and 

Fs  =  thermal  resistance  of  the  film. 


[W/m2] 

:  C] 

°C] 

;w/m°C] 

>] 

W/m°C] 

;j/kg°C] 

k^/m3] 


This  expression  is  valid  only  when  the  heat  is  conducted  through  the  film  and 
convected  away  by  the  discs,  a  condition  that  is  true  for  most  concentrated 


contacts.  The  strength  q  of  the  source  is  given  by: 
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(3-11)  q  =  rAU  . 


The  thermal  resistance  of  the  film  Fs  may  be  found  by  calculating  the  average 
film  temperatures  due  to  distributed  sources  and  is  normally  taken  to  be: 


(3-12)  Fs  =  0.1. 


In  order  to  proceed  any  further  we  need  a  relationship  between  temperature 
and  the  shear  strength  of  the  fluid.  A  typical  relationship  that  has  its 
roots  in  the  Eyring  theory  of  fluid  transport  is  given  by: 


(3-13) 


t{8)  = 


1 

Z 


A  +  BP  +  (0+D)  Ln 


where, 


A  =  viscosity  temperature  constant,  [°C] 

B  =  pressure  viscosity  constant,  [°C/Pa] 

C  =  non-linear  shear  stress  constant,  [°C/Pa] 

D  =  fluid  solidification  temperature,  [°C] 

E  =  fluid  viscosity,  and  [Pa.s] 

P  =  pressure  of  the  fluid  in  the  contact.  [Pa] 


At  first  sight  it  seems  that  this  equation  has  five  disposable  constants  in 
it.  However,  two  of  these  constants  (A  and  D)  may  be  obtained  from  the 
atmospheric  viscosity  temperature  relationship,  and  one  more  (B)  can  be 
obtained  from  the  Barus  viscosity  pressure  relationship.  The  constant  D  is 
known  as  the  solidification  temperature,  the  temperature  to  which  the  fluid 
should  be  cooled  to  become  solid  like  under  atmospheric  pressure.  Only  the 
constants  C  and  E  need  to  be  determined  experimentally  from  the  traction 
results  and  this  will  be  done  in  the  next  chapter.  It  should  be  noted  here 
that  the  ultimate  aim  is  to  derive  the  fluid  traction  parameters  such  that 
they  apply  for  all  the  experimental  conditions  reported  here. 

By  using  equations  (3-9),  (3-10),  (3-11),  and  (3-13)  the  average  thermal 
shear  stress  can  be  obtained  for  a  given  set  of  conditions.  In  equation  (3-8) 
we  need  the  ratio  of  the  average  contact  shear  stress  under  thermal  conditions 
to  that  under  isothermal  conditions.  This  is  really  the  ratio  of  the  shear 
stress  given  by  equation  (3-13)  evaluated  at  the  shearplane  temperature  6C 
(from  equation  3-9)  and  the  shear  stress  as  evaluated  at  the  inlet  temperature 
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conditions  60.  The  other  contact  conditions  remain  the  same  for  this  ratio 
calculation. 

The  method  that  is  outlined  above  was  used  for  the  traction  data  analysis 
as  reported  by  Tevaarwerk  [1985].  For  the  analysis  of  the  traction  data  as 
reported  here,  several  improvements  were  made  in  the  model  and  these  changes 
are  outlined  in  the  next  chapter. 
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4-0  EXTRACTION  OF  THE  TRACTION  PARAMETERS 

In  order  to  use  the  traction  model  to  calculate  the  level  of  traction  in 
a  given  contact  condition,  the  relevant  traction  parameters  need  to  be 
extracted.  In  equation  (3-1)  we  saw  that  both  the  elastic  effects  and  the 
nonlinear  viscous  effect  account  for  the  behavior  of  the  fluid.  Elastic 
effects  are  typically  felt  under  conditions  of  high  contact  pressure,  low 
temperature,  and  high  speed. 

4-1  Extraction  of  The  Shear  Modulus  Parameters 

As  discussed  in  chapter  3  the  initial  linear  slope  on  a  traction  curve 
can  be  the  result  of  either  a  viscous  or  an  elastic  response  of  the  material 
in  the  contact  to  the  strain  implied.  The  parameter  that  determines  the 
actual  response  is  the  dimensionless  grouping  of  the  relaxation  time  of  the 
material  in  the  contact  and  the  transit  time  of  this  material  through  the 
contact.  For  a  simple  Maxwell  type  material  this  number  is  known  as  the 
Deborah  number  and  is  given  by: 


(4-1)  De  -  , 


where, 


17  =  viscosity  of  the  fluid  in  the  contact, 

U  =  transit  velocity  of  fluid  in  the  contact, 
a  =  semi  contact  length  in  rolling  direction,  and 
G  =  shear  modulus  of  the  fluid  in  the  contact. 


[Pa.s] 

[m/s] 


Under  the  assumptions  of  a  simple  pressure  distribution  according  to  Hertz,  a 
constant  film  thickness  in  the  contact,  and  a  constant  shear  modulus  over  the 
contact  area,  the  initial  small  strain  behavior  for  the  material  will  be 
elastic  if  the  Deborah  number  is  larger  than  10  and  will  be  viscous  in 
response  if  the  Deborah  number  is  less  than  0.1.  In  between  these  two  values 
the  response  is  due  to  a  mixed  viscoelastic  behavior. 

Many  assumptions  have  been  made  before  arriving  at  this  point,  however 
even  with  more  complicated  analyses  where  the  pressure  is  allowed  to  influence 
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the  viscous  and  elastic  properties  it  is  found  that  the  transition  points 
occur  at  about  the  same  Deborah  numbers.  Also  for  the  range  of  parameters 
normally  encountered  in  traction  contacts  the  Deborah  number  is  such  that  the 
initial  slope  of  the  traction  curve  is  almost  always  governed  by  the  elastic 
properties  of  the  contacting  material.  In  the  analysis  of  the  data  as 
performed  here  this  assumption  is  implicit. 

4-1-1  Shear  Modulus  For  Constant  Properties 

When  the  initial  linear  response  is  completely  elastic  it  is  quite  easy 
to  calculate  the  value  for  the  actual  modulus  that  caused  this  slope.  Under 
the  assumptions  of  constant  properties  throughout  the  contact  this  modulus  can 
be  extracted  from  the  initial  slope  using  the  following  equation: 


(4-2)  G  =  irm  — ^  [Pa]  . 

This  equation  is  applicable  regardless  of  whether  the  slope  is  measured  under 
longitudinal  slip  or  under  side  slip.  Using  the  known  values  of  m,  P0,  and  a, 
and  a  calculated  value  of  h  for  film  thickness,  it  is  quite  simple  to 
calculate  the  fluid  shear  modulus  G. 

4-1-2  Shear  Modulus  With  a  Simple  Compliance  Correction 

One  of  the  criticisms  that  is  often  raised  at  the  above  analysis  is  that 
it  neglects  the  influence  of  the  disk  compliance  on  the  measured  slope.  Disk 
compliance  is  the  result  of  the  elastic  creepage  of  the  disk  material  due  to 
the  tractive  stresses  on  the  surface.  The  traction  response  in  the  initial 
linear  range  is  affected  by  this  disk  creepage  in  that  it  makes  the  slopes 
lower  than  if  the  discs  were  infinitely  stiff.  An  exact  correction  of  the 
modulus  for  the  disc  compliance  is  not  possible  at  the  moment.  The  analysis 
that  is  presented  here  is  that  due  to  Johnson  and  Roberts  [1974]. 

If  we  let  m‘  be  the  slope  of  the  traction  curve  for  dry  rolling  bodies, 
then  from  the  addition  of  the  compliances  of  the  discs  and  the  film  a  simple 
corrective  term  for  the  shear  modulus  may  be  derived  as  shown  in 
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Equation  (4-3).  From  this  expression  it  is  obvious  that  when  the  measured 
slope  approaches  the  dry  slope  the  corrected  value  for  the  shear  modulus  tends 
to  infinity. 


(4-3)  ^  =  6(5^  ,  M 

where  Gc  =  simple  compliance  corrected  modulus. 

The  dry  slope  m'  can  be  calculated  from  the  expression  given  by  Kalker  [1967] 
as: 


(4-4)  m1  -  ,  [-] 

where, 

Gs  =  shear  modulus  of  the  disc  material,  and  [Pa] 

Q  =  Kalker  coefficient.  [-] 

The  value  of  the  Kalker  coefficient  depends  on  the  aspect  ratio  of  the 
contact  and  the  direction  ot  slip.  Below  is  a  table  which  gives  these  values 
for  the  tests  as  reported  here. 


Table  4-1.  Kalker  Coefficients  Q  for  Sideslip 

aspect  ratio 

coefficient 

k 

Q 

1 

.56 

2 

.70 

5 

.81 

23 


4-2  Extraction  of  the  Large  Strain  Parameters 

The  large  slip  region  ,  that  is  the  region  beyond  the  traction  peak,  is 
almost  exclusively  governed  by  the  dissipative  element  in  the  rheological 
equation.  With  the  elastic  effects  becoming  insignificant  it  is  now  quite 
easy  to  extract  the  governing  parameters  for  this  region.  In  essence  what  is 
required  is  a  reverse  analysis  of  the  traction  calculation  normally  used  for 
the  calculation  of  the  traction  curves.  When  the  elastic  response  of  the 
fluid  is  no  longer  dominant,  the  shear  strain  rate  in  the  film  is  taken  by  the 
viscous  element  and  we  may  express  equation  (3-1)  as: 

(4-5)  F (r)  =  7  =  ^  sinh  (  ) 


When  the  argument  r/r,  is  larger  than  1  this  may  be  further  simplified  to 
the  following: 


(4-6)  F(t)  =  7  =  ^  Exp  (  -^  )  . 

This  expression  relates  the  stress  to  strain  response  for  a  small  element 

O 

of  fluid  that  is  subjected  to  a  shear  strain  rate,  7,  at  a  given  pressure  and 
a  given  temperature.  Each  of  the  parameters  in  this  equation  can  be  a 
function  of  temperature  or  pressure  or  both.  In  the  analysis  of  traction  it 
is  normally  sufficient  to  treat  the  fluid  viscosity  as  a  function  of 
temperature  and  pressure  and  to  make  the  non-linear  stress  parameter  a 
function  of  temperature  only.  The  following  relationships  will  be  used  here: 

---  for  the  viscosity  we  will  combine  the  Vogels  temperature  viscosity  and 
the  Barus  pressure  viscosity  relationships  as  follows: 

r,  ( 6 ,  p)  =  %  exp  +^) 


(4-7) 


[Pa.s] 


24 


where, 


Vo 

=  viscosity  constant, 

[Pa.s] 

A 

=  temperature  viscosity  constant, 

[8C] 

Dy 

=  solidification  constant  for  Vogels  equation, 

[•c] 

B 

=  pressure  viscosity  constant, 

[°C/Pa] 

Dp 

=  solidification  constant  for  Barus  equation, 

[°C] 

6 P 

=  local  temperature,  and 

[°C] 

P 

=  local  pressure. 

[Pa] 

—  for  the  temperature  dependence  of  rs  on  the  local  temperature  we  will 
use  the  following: 


(4-8) 

where, 


Q  +  D. 


I 


C  =  non  linear  shear  stress  constant  [°C/Pa] ,  and 
Ds  =  solidification  constant  for  Eyring  equation  [°C] . 


These  parameters  (except  for  the  temperature  and  the  pressure)  are  all 
constants  for  a  given  fluid. 

In  order  to  proceed  any  further  we  have  to  make  some  assumptions  about 
the  conditions  within  the  contact.  First,  we  will  assume  that  average 
conditions  exist  for  the  contact.  Hence  temperature  and  pressure  are  only 
mean  temperatures  and  pressures.  Second,  we  will  assume  that  the  important 
temperature  for  the  analysis  is  the  shear  plane  temperature  as  given  by 
Equation  (3-9).  Third,  the  shear  strain  rate  7  will  be  assumed  to  be  constant 
throughout  the  film  and,  in  the  case  of  the  traction  experiments,  is  constant 
over  the  contact  area.  Its  magnitude  is  given  by: 


(4-9)  7  =  Av/h,  and  [s'1  ] 

h  =  central  f i lmthickness  [m] 

The  central  film  thickness  is  calculated  from  the  expressions  by  Hamrock 
and  Dowson  [1977].  Further  modifications  to  this  film  thickness  to  allow  for 
inlet  shear  heating  are  made  by  using  the  Murch  and  Wilson  [1975]  approach. 
Hence  in  total  we  will  make  the  following  simplifications: 
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(4-10) 


local 


average  equation  # 


temperature 

e 

6C 

(3-9) 

shear  strain  rate 

0 

7 

Av/h 

(4-9) 

pressure 

P 

2P0/3 

(2-1) 

shear  stress 

r 

Fy/*ab 

By  substituting  the  expressions  from  equation  (4-7)  to  (4-9)  into  Equation 
(4-2)  we  obtain: 


(4-11) 


Av  ec  +  Ds 

TT=  — ZCV7 


_  (  Tc  c  A  B  P  \ 

Exp  \“57R)5  ■  '  “gr+Dp  v 


At  a  first  glance  there  appear  to  be  7  parameters  that  can  be  used  to  fit  the 
data.  This  is  not  the  case  however  because  several  of  these  are  determined  by 
the  results  from  other  experiments.  For  example,  the  correlation  of  the 
viscosity  data  for  the  lubricants  to  equation  (4-3)  solves  for  the  constants 
A,  B,  Dv,  and  Dp.  It  remains  therefore  to  find  the  constants  Ds  and  C.  These 
can  be  derived  from  the  thermal  region  of  the  traction  curve  itself  by  curve 
fitting  Equation  (4-11)  to  it.  For  this  purpose  it  is  better  to  write  this 
equation  in  a  slightly  different  form: 


(4-12) 


BP 


+  Ln 


( 


2CV„  av 

M.)  h 


From  the  above  equation  it  is  apparent  that  if  this  relationship  holds  then 
the  results  from  the  traction  measurements  should  form  a  straight  line  when 
plotted  in  the  above  fashion.  The  slope  of  this  line  gives  the  value  for  C 
while  the  value  of  Ds  can  be  calculated  by  adjusting  it  until  the  best  fit  is 
obtained. 
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4-3  Reduced  Pressure  Effects 


When  a  fluid  film  is  present  in  the  contact  zone  the  pressure 
distribution  is  not  strictly  Hertzian  but  modified  by  the  fluid.  Due  to  the 
hydrodynamic  action  the  pressure  distribution  is  more  peaky,  and  is  spread 
over  a  broader  area.  This  will  lead  to  a  reduction  in  the  mean  pressure  in  the 
contact.  The  reduced  pressure  may  be  calculated  from  the  pressure  ratio  as 
given  below: 

(4-13)  Pr  =  P  [1  -  4  Ge-0  25  YI'3  e('2  3/k)]  , 

where, 

Ge  =  Johnson’s  elasticity  constant,  [-] 

P  =  mean  Hertz  contact  pressure,  [Pa] 

YI  =  Inlet  shear  heating  factor,  [-] 

k  =  aspect  ratio  (b/a),  [-] 

Pr  =  reduced  contact  pressure,  [Pa] 


4-4  Transit  Time  Dependence 

The  fluid  in  the  contact  is  subjected  to  very  high  pressure  gradients  for 
a  very  short  time  duration.  It  is  therefore  not  likely  that  the  viscosity  in 
the  contact  has  reached  equilibrium  conditions.  How  near  equilibrium  that  it 
is  will  depend  on  several  factors  such  as  the  time  it  is  subjected  to  the 
pressure  and  the  ease  with  which  the  molecules  can  rearrange  themselves. 

A  suitable  dimensionless  grouping  that  takes  these  various  factors  into 
account  is  given  as: 


(4-14)  *  =  , 

where, 

$  =  dimensionless  time  delay  parameter,  [-] 


U  =  mean  roller  speed,  [m/s] 

J?(0O)  =  viscosity  at  inlet  temperature,  [Pa - s] 

a  =  semi  Hertz  contact  length,  and  [m] 

P0  =  Hertz  contact  pressure.  [Pa] 


27 


This  parameter  is  used  in  conjunction  with  the  viscosity  constant  in  the 
following  way: 


(4-15)  V0  =  Efi*2,  [Pa.s] 


where, 

E0  =  equilibrium  viscosity,  and  [Pa.s] 

E2  =  time  delay  exponent.  [-] 


The  single  constant  V0  has  now  been  replaced  by  this  two  constants  equation 
so  we  have  increased  our  degree  of  freedom  by  one  for  the  system. 

4-5  Thermal  Conductivity  Temperature  Dependence 

It  may  be  expected  that  the  thermal  conductivity  of  the  fluid  is  a  function  of 
the  temperature.  This  has  a  direct  influence  on  the  thermal  resistance  of  the 
fluid  in  the  contact.  In  order  to  take  this  effect  into  account  the  thermal 
conductivity  is  fitted  to  the  following  equation: 

(4-16)  kf  =  '  [W/m  °C] 

Since  the  shear  plane  temperature  occurs  in  equation  (4-12)  and  the  thermal 
conductivity  occurs  in  equation  (3-9),  the  solution  to  both  the  thermal 
conductivity  and  the  shear  plane  temperature  must  be  done  on  an  iterative 
basis. 

4-6  Determination  of  the  Experimental  Parameters 

We  are  now  in  a  position  to  use  the  data  obtained  from  the  traction 
experiments  and  fit  the  rheological  model  to  them.  As  stated  earlier,  the 
viscosity  temperature-pressure  data  for  the  traction  data  reduction  may  be 
obtained  from  isobaric  and  isothermal  viscosity  measurements.  Similarly  the 
thermal  conductivity  data  can  be  obtained  from  published  information  on  RP1  in 
the  Liquid  Propellant  Manual  [1982].  The  following  sections  will  deal  with 
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the  determination  of  the  actual  numerical  values  of  the  various  constants 
needed  for  the  equations  of  state  and  the  constants  needed  for  the  fluid 
rheology. 


4-6-1  Determination  of  Viscosity  Temperature  Pressure  and 
Temperature  Conductivity  Parameters 

Not  all  of  the  constants  in  Equation  (4-12)  need  to  be  determined  from 
the  experimental  traction  data.  The  temperature  viscosity  characterization 
may  be  used  to  solve  for  A  and  Dv  while  the  pressure  viscosity  effects  can  be 
used  to  solve  for  B  and  Dp. 

The  viscosity  temperature  data  may  be  obtained  from  the  Liquid  Propellant 
Manual  [1982].  Pressure  viscosity  data  for  RP1  was  taken  from  the  data 
generated  by  Bridgeman  [1970]  for  kerosene.  Using  these  two  sources  of 


information  allows  the  determination 

of  the  constants 

in  Equation  (4-7) 

(4-7) 

V  (9.  p)  =  Vo  exP 

(  A  ,  pB  ) 

V  FIT  +  6RT  / 

,  [Pa.s] 

where, 

v  P 

A  =  542.4, 

[°C]  , 

B  =  544, 

[GPa-1] 

Dy  =  143, 

[°C] 

D  =  25.6,  and 

[°C] 

V0  =  6.33  x  10’6 

[Pa-s] 

Thermal 

conductivity  data  for  RP1  was 

taken  from  the  Vargaftik  [1983]. 

Fitting 

this  data  to  Equation  (4-16), 

(4-16) 

k' 

kf  dc  +  k"  * 

[W/m  °C] 

results 

in  the  following  values  for  k 

'  and  k" 

k'  =  61.6,  and 
k"  =  535.5 


[W/m  °C] 
[°C] 
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4-6-2  Experimental  Shear  Modulus  from  Traction  Data 

Similar  to  the  non-linear  stress  parameter  rc<  it  is  found  that  the  fluid 
shear  modulus  as  extracted  from  traction  tests  is  influenced  by  the  inlet 
temperature,  by  the  contact  pressure,  and  by  the  rolling  speed.  As  a  fairly 
simple  relationship  the  following  expression  is  used: 


(4-17) 

where, 


Gj  +  G2  P  +  G3  |  U 
0  +  Dy 


Gj  3  =  experimental  constants, 

P  '  =  contact  pressure, 

U  =  rolling  velocity,  and 

Dv  =  Vogels  fluid  solidification  temperature. 


[Gpa] 


[GPa] 

[m/s] 

[°C] 


Using  the  initial  slopes  from  the  traction  experiments  and  some  regression 
analysis  the  following  values  were  found  for  the  constants  Gj  3. 


G,  1.0 

G2  =  0.25 

Gg  =  -0.005 


[GPa  »C] 

[°C] 

[GPa  °C  s/m] 


To  obtain  this  data,  the  high  pressure,  low  temperature  curves  were  mostly 
used  in  the  calculations. 


4-6-3  Experimental  Traction  Constants 


By  combining  the  Equations  (4-8),  (4-9)  and  (4-10)  the  large  slip  traction 
data  may  be  fitted  to  the  following  equation: 


(4-12) 


1 

Z  L 


BP 


+  Ln 


( 


2CV„  Av 

M.)  !>' 


I 


where  the  shear  plane  temperature  is  calculated  from  equation  (3-9)  as: 
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(3-9)  ec  -  60  =  q  [f,  + 

kf  (  JPe 

P8  is  given  by  Equation  (3-10)  as: 


The  constants  that  remain  to  be  solved  for  at  this  point  are: 

C(  Ds,  E0,  and  E2. 

These  can  now  be  solved  for  by  using  a  successive  regression  analysis  on  the 
large  slip  traction  data.  This  was  performed  using  the  experimental  traction 
data  as  reported  in  Chapter  2.  A  summary  of  the  test  data  used  and  the 
calculated  values  for  some  of  the  intermediate  parameters  are  shown  in 
Table  4-2.  The  thermal  properties  for  the  test  disc  material  was  taken  to  be: 

ks  =  15,  [W/m  CC] 

(4-18)  Cs  =  500.  [J/kg  °C] 

The  degree  of  fit  obtained  from  the  data  on  Equation  (4-12)  is  shown  in  Figure 
4-1.  From  this  regression  analysis  the  following  constants  resulted: 

C  =  2.75  x  10'5-  [°C/Pa] 

Ds  =143,  [°C] 

E  =  6.96  x  10'5,  and  [Pa.s] 

E2  =  -0.057.  [-] 


(4-19) 
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TABLE  4-2.  TRACTION  CURVE  DATA  ANALYSIS  FOR  RP1 


The  following  oil  properties  were  used  in  the  analysis- 

Viscosity  temperature;  A=  542.4  ’C  D=143.0  ’C  V0=  .0000B33  Pas 

B=  544.0  1 C/GPa  Dp  =  25. B  ’C 
K  ’  •  535.5  U/m’C  K”=  61.7  ’C 
Peak  Factor=  .98 


Viscosity  Pressure 
Thermal  Conductivity 
Eyring  Dpress=  25.6 


Disc  Thermal  Conductivity  =  Ks=  15  U/m’C  Specific  Heat=500  J/kg’C 


’*»#4*t*******«t«*>ft*«ft*»*«#»*»«^ 


►#♦*#######*#***#*#♦***♦ 


TEST  INCLUDED  IN  THE  MULTIPLE  TRACTION  DATA  ANALYSIS 


Test  # 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

T  shear 

H0 

Yi 

Plot 

( - ) 

( GPa  ) 

( m/  s  ) 

(  ’  C  > 

(N  ) 

( -  > 

(  mm  ) 

(  -  > 

(  ’C  ) 

mi  c 

— 

Mar  k 

RP196M1TR 

1 

.01 

10 

.3 

34 

200 

1.7 

.235 

n 

L. 

47 

.097 

.98 

A 

RP196M2TR 

1 

.01 

30 

.4 

35 

200 

1.7 

.235 

2G 

61 

.193 

.92 

B 

RP196M3TR 

1 

.01 

49 

.6 

39 

200 

1.7 

.235 

65 

71 

.237 

.87 

C 

RP196M4TR 

1 

.27 

11 

.5 

37 

400 

1.7 

.297 

n 

C 

61 

.093 

.98 

D 

RP196M5TR 

1 

.27 

30 

.8 

38 

400 

1.7 

.297 

18 

76 

.  175 

.92 

E 

RP19GMGTR 

1 

.27 

50 

.0 

41 

400 

1.7 

.297 

44 

92 

.211 

.87 

F 

RP19SM7TR 

1 

.60 

11 

.9 

39 

800 

1.7 

.374 

1 

75 

.086 

.98 

G 

RP196M8TR 

1 

.60 

30 

.5 

40 

800 

1.7 

.374 

11 

99 

.  150 

.93 

H 

RF19GM9TR 

1 

.60 

50 

.5 

44 

800 

1.7 

.374 

31 

117 

.  188 

.88 

I 

RF19GMATR 

1 

.92 

11 

.8 

40 

1380 

1.7 

.448 

1 

91 

.077 

.98 

J 

RP196MBTR 

1 

.92 

31 

.4 

44 

1380 

1.7 

.448 

7 

119 

.  128 

.94 

K 

RP19GMCTR 

1 

.92 

50 

.0 

49 

1380 

1.7 

.448 

19 

141 

.157 

.89 

L 

RP19GN1TR 

1 

.01 

49 

.5 

58 

200 

1  .7 

.235 

37 

84 

.175 

.90 

M 

RP19GN2TR 

1 

.01 

30 

.7 

57 

200 

1.7 

.235 

13 

74 

.135 

.94 

N 

RP196N5TR 

1 

.27 

30 

.  9 

60 

400 

1.7 

.297 

9 

87 

.123 

.94 
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FIGURE  4-1.  DEGREE  OF  FIT  OF  RP1  TRACTION  DATA  TO  EQUATION  (4-12) 
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5-0  TRACTION  PREDICTION 

The  ultimate  aim  in  obtaining  the  fluid  rheological  properties  is  to 
solve  for  the  traction  response  of  a  concentrated  contact  under  general 
conditions  of  slip,  spin,  and  geometry.  In  this  report  we  will  restrict 
ourselves  to  the  conditions  of  sideslip  or  longitudal  slip  only.  Calculation 
methods  for  other  contact  kinematic  conditions  are  outlined  in  Tevaarwerk 
[1979a]  and  Tevaarwerk  [1979b]. 

As  another  way  to  observe  the  faithfulness  of  the  traction  modeling  it 
may  be  instructive  to  predict  the  measured  traction  curves  on  the  reduced  data 
as  extracted. 

5-1  Comparison  Between  Measured  and  Predicted  Traction 

The  process  essentially  involves  the  reverse  operation  to  that  used  for 
the  extraction  of  the  traction  parameters.  Starting  with  Equation  (3-8) 


(3-8)  J4t  =  J4e  +  J4p 


where,  J4e  = 


4  S 

ir(l+S2)2 


2 

r 


tan_1S  + 


S(S2  -1) 
(1+S2)2 


,  and 


s  =  2  J1  =  2  G  a  AU 

3  Ik  3  r0  h  U 

As  an  average  shear  modulus  G  we  have  to  use  the  combination  of  fluid 
shear  modulus  and  the  elastic  compliance  of  the  discs.  This  apparent  modulus 
can  be  calculated  from  the  compliance  corrected  modulus  as  given  in 
Equation  (4-17)  together  with  the  relationships  (4-2),  (4-3),  and  (4-4).  This 
results  in: 
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(5-1) 


G  = 


1  + 


Gc 

8Qa  Gc 
3rhGs 


[Pa] 


where, 

Gc  =  compliance  corrected  modulus  as  given  by  Equation  (4-17), 

a  =  semi  contact  size  in  rolling  direction,  [m] 

Gs  =  shear  modulus  of  disc  material,  [Pa] 

Q  =  Kalker  coefficient,  and  [-] 

h  =  central  contact  film  thickness.  [m] 


The  average  fluid  shear  strength  r0  based  on  the  inlet  temperature  can  be 
calculated  directly  from  Equation  (4-12)  by  using  90  as  the  shear  plane 
temperature, 


(5-2) 


BP 


+  Ln 


( 


2CV.  a v 

M.)  h' 


) 


,  and 


using  the  fluid  constants  as  given  in  Equation  (4-19). 

To  solve  for  the  ratio  of  fluid  shear  strength  at  the  shear  plane  temperature 
6C  and  the  inlet  temperature  90  requires  a  simple  iteration  since  0C  is  a 
function  of  rc  and  vice-versa  as  indicated  by  Equations  (3-9),  (3-10),  (3-11) 
and  (4-12).  Additional  to  the  temperature  dependence  of  rc  through  dc  there 
is  also  the  influence  of  thermal  conductivity.  In  most  cases  only  two  or 
three  iterations  are  needed  for  convergence.  By  using  this  technique,  and  the 
rheological  data  as  extracted,  the  sideslip  traction  can  be  predicted  for  each 
of  the  measured  traction  curves.  The  comparison  between  measured  and 
predicted  traction  may  be  observed  from  the  graphs  included  in  Appendix  II. 
Overall  the  degree  of  predictions  correlates  very  well  with  the  experimentally 
observed  results  indicating  a  high  degree  of  confidence  in  the  traction 
model ing. 
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5-2  Traction  Predictions  for  Line  Contact 

Traction  prediction  for  line  contact  is  essentially  the  same  as  that 
outlined  in  Section  5-1.  The  fundamental  fluid  parameters  are  independent  of 
the  exact  contact  geometry  so  there  would  be  no  change  there.  Changes  do  need 
to  be  made  in  the  following: 

1)  Film  thickness  calculations, 

2)  Kalker  correction  factor, 

3)  Mean  contact  pressure,  and 

4)  Shear  plane  temperature  calculation. 

The  film  thickness  calculations  should  reflect  the  increased  film  that  is 
found  when  rollers  form  a  line  contact.  In  the  line  contact  calculations 
performed  here,  the  central  film  thickness  calculations  according  to  Dowson- 
Higginson  [1977]  are  used.  The  effect  of  the  increased  film  thickness  is  a 
slight  reduction  in  the  traction  in  the  thermal  region  because  of  increased 
thermal  resistance  of  the  thicker  film  as  compared  with  elliptical  contacts. 

Kalker's  correction  factor  for  the  compliance  of  the  discs  should  now 
reflect  the  direction  of  the  slip  and  the  contact  geometry.  From  Kalker 
[1967]  we  find  that  a  value  of  Q  for  line  contact  is: 

(5-3)  Q  =  0.6.  [-] 

When  the  traction  data  was  fitted  to  the  fluid  rheological  model,  averaged 
quantities  were  used.  The  average  pressure  in  a  line  contact  is  somewhat 
higher  for  a  given  peak  contact  pressure  than  in  an  elliptical  contact.  From 
Hertz  stress  considerations  we  can  calculate  a  correction  factor  for  line 
contact.  Hence  for  the  same  peak  Hertz  pressure  we  should  use  a  mean  pressure 
that  is  about  18  percent  higher  in  the  traction  calculations. 

The  shear  plane  calculations  are  based  on  Equation  (3-9).  This  equation 
consists  of  two  thermal  resistance  terms,  1)  that  due  to  the  film,  and  2)  that 
due  to  the  contact.  The  change  in  film  resistance  is  taken  care  of  by  using 
the  film  thickness  calculations  for  line  contact.  The  change  in  contact 
resistance  can  be  corrected  for  by  using  the  ratio  of  contact  resistance  for 
circular  contact  and  rectangular  contact.  Jaeger  [1942]  gives  this  ratio  as: 


This  is  based  on  average  contact  temperature.  Using  the  corrected  terms  in 
relevant  equations,  calculations  identical  to  these  for  elliptical  contact 
under  side  slip  can  be  carried  out. 

5-2-1  Line  Contact  Traction  Predictions 

One  of  the  aims  of  this  investigation  was  to  provide  traction  data  on  RP1  that 
can  be  used  in  the  bearing  analysis  code  SHABERTH.  The  traction  data  required 
by  SHABERTH  is  in  the  form  of  a  number  of  selected  points  from  line  contact 
traction  test.  The  above  outlined  technique  was  used  to  modify  the  data 
obtained  using  elliptical  contact  traction  tests  into  line  contact  data. 

Using  this  method  the  following  line  contact  traction  test  curves  were 
predicted. 


Contact  Pressure 

P0  = 

1.6,  2.3,  3.0, 

[GPa] 

Rolling  Speed 

U  : 

30,  60,  100, 

'm/s] 

Temperature 

5  o' 

40,  75,  100,  and 

[°c] 

Slide/Roll  Ratio 

AU. 

“IT- 

0.0  to  0.40,  | 

[-] 

The  predicted  line  contact  traction  curves  are  included  in  Appendix  III. 
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APPENDIX  I 


ORIGINAL  MEASURED  TRACTION  DATA 
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APPENDIX  II 


COMPARISON  BETWEEN  MEASURED  AND 
PREDICTED  ELLIPTICAL  CONTACT  DATA 
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APPENDIX  III 

PREDICTED  LINE  CONTACT  DATA 
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